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P. Donohue^ and T. Giamarchi* 

Laboratoire de Physique des Solides, CNRS-UMR 85002, 
Universite Paris-Sud, Bat. 510, 91405 Orsay, France 



o 
o 

(N 
S3 

a 
in 



(Dated: February 6, 2008) 

We show that in a commensurate bosonic ladder, a quantum phase transition occurs between a 
Mott insulator and a superfluid when interchain hopping increases. We analyse the properties 
of this transition as well as the physical properties of the two phases. We discuss the physical 
consequences for experimental systems such as Josephson Junction arrays. 
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The Mott transition is one of the most striking conse- 
quences of interactions in a quantum system. Although 
mostly studied in fermionic systems, it is even more spec- 
tacular in interacting bosonic ones where the transition 
occurs between a Mott insulating state where bosons are 
well localized on the lattice and a superfluid state, with 
zero resistivity. The existence of such a Mott transi- 
tion for bosons has been proven in one dimensioncl us- 
ing bosonization and studied in higher dimensions us- 
ing a scaling theorycl. Since then the physical properties 
of this transition have been intensely investigated^. For 
this problem, one dimension is a specially interesting case 
since: (i) interactions play an extremely important role 
leading, even for bosons, to a Luttinger liquid state; (ii) 
no true superfluid condensate exists. Nevertheless a true 
Mott-superfluid transition takes place whose critical and 
physical properties are now well understood both at potft- 
mensurate densities and close to commensurabilitjU'QQ. 
The existence and characteristics of this transition have 
been investigated by various analytical and numerical 
techniqueaj'ErLlQ an d experimentally in Josephson Junc- 
tion (JJ) arraysBEl 

Even if the properties of a single bosonic chain are 
now rather well understood, a crucial question is how 
the characteristics of the transition evolve when going 
from one to two dimensions. Indeed if one considers cou- 
pled bosonic chains there is clearly a competition between 
the interactions (thus the Mott phase) and the interchain 
hopping that tends to delocalize the system and push it 
to the superfluid phase. One way to tackle these ques- 
tions analytically is to focus on a small number of coupled 
chains, i.e. on ladders. Such ladder systems have also ex- 
perimental realizations in JJ arrays. A similar approach 
was followad.in.ici'mionic systems for which related ques- 
tions aroseOll3ll3. For commensurate fermionic ladders, 
quite surprisingly, only a strong crossover is found but 
the system remains insulating regardless of the strength 
of the interchain hopping. This prompts for the question 
of the existence of the Mott-superfluid transition in the 
bosonic ladder. So far the physical properties of bosonic 



ladder|S-ihave been studied only for the incommensurate 
laddeili-3 jar for a commensurability of one boson every 
two sitesEj for which the ladder-is essentially equivalent 
to an anisotropic spin Iaddenl3't3. However the case of 
one boson per site, where the commensurability being 
higher leads to a stronger and more involved competi- 
tion between the Mott phase and the interchain hopping 
is still to be understood. 

We thus investigate in the present paper a bosonic lad- 
der, at or close to the commensurate filling of one boson 
per site. We show that for this system there is a true 
transition, when interchain hopping is increased, between 
a Mott insulator and a superfluid phase. Large coupling 
expansions show that this transition is in the Beresinskii- 
Kosterlitz-Thouless (BKT) universality class at commen- 
surate filling and in the commensurate-incommensurate 
one for small doping. Various correlation functions such 
as the single particle one, show universal power law de- 
cay at the transition. The transport properties, that we 
compute can provide a check of the existence of this tran- 
sition in experimental systems such as the JJ arrays. 

We start with the Hamiltonian 



H = -t yfoj, a &i+i,q + h-c.) + — y^n^ a (rtf |Ct - 1) 
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where a = 1, 2 is the chain index, t and t± are respec- 
tively the intra and interchain hopping, U is the on-site 
particle interaction, and p is the chemical potential. In 
the following, unless we specify otherwise, p is chosen to 
impose one boson per site. To describe the low energy 
properties of (jl|) it is convenient to use the density-phase 
representation of the bosonslil: 

p{x) w po - -d x <f> + po cos(2</>) (2) 
A^^oce^)^) ( 3 ) 
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-d x <Ky),9(y') 



= iS(y' - y), 
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where a is a short distance cut-off of the order of the 
lattice spacing. Using these variables the single chain 
Hamiltonian can be rewritten as: 



H = H° 



9u 



(27ra) 



H°(u,K) = — dx 



dx cos 2(f> 
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All interaction effects are hidden in the Luttinger liquid 
parameters u, the sound velocity, and K a coefficient that 
controls the asymptotic decay of the correlation func- 
tions. Free bosons correspond to K — oo and hard core 
ones (U = oo) to K = 1. Note that for one boson per 
site, the U — oo case is a trivial band insulator. The 
cos(2(/>) term where g u oc U describes the umklapp scat- 
tering of the bosons,on the lattice and is responsible for 
the Mott transitionQ Note that the continuous form (|^) 
is much more general than the microscopic hamilonian 
([j]) and describes systems with longer range interactions 
as well. Using the symmetric basis S)O = {cf>\ ± 02)/v2) 
(@) gives: 

H = H° a {u,K a )+H° s (u,K s ) 



(27ra) 2 j 

■ — / dx cos(V20 a 
na J 



dx cos 
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where H® a are denned by (^). For (|l|) one has u s . a — u 
and A*. = K. 

In (ffi) there is a competition between the umklapp 
term, that favors localization of the charge hence order 
in the field 4> a (and <p s ) and the interchain coupling that 
wants to order the relative superfluid phase 9 a between 
the two chains. Since these two fields are conjugate one 
can naively expect a transition between the two types of 
order for the antisymmetric field. The resulting compe- 
tition affects the symmetric modes, inducing a metal- 
insulator transition. To investigate this transition we 
use a renormalization group (RG) method in powers of 
the umklapp and interchain hopping. The RG equations 
read: 



(2 - KJ2 - K s /2)g u 



-± = (2-K- 1 /2)t ± 



dl 
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where g u — g u /u and t± = (4ir 2 a)t±/u are dimension- 
less coupling constants. For 1/4 < K < 2, (^) show that 
both t± and g u are relevant operators that thus tend 



to order 9 a and (f> a respectively. Two different behav- 
iors occur depending on which operator becomes of or- 
der unity first: (i) if the umklapp scattering eventually 
dominates the flow it leads to ordering in <f) s and <p a and 
we recover an insulating behavior similar to the one of 
a single chain; (ii) if the interchain hopping dominates 
then 9 a orders and the umklapp scattering becomes ir- 
relevant. The only term that can localize the symetric 
component is now a term generated to second order of 
the form cos(v / 80 s )- This operator has a dimension 2K 
instead of K for a single chain umklapp, and is thus much 
less relevant. The precise interaction for which this op- 
erator is relevant depends on its amplitude that can only 
be computed reliably when g„ < 1. Even if this is not 
the case, one can show in the large t± expansion below 
that there is a range of interactions between t/ s c inglo chain 
and U = oo for which this operator is indeed irrelevant. 
Thus at fixed interactions a transition between the Mott 
insulator and a superfluid for a given value of the in- 
terchain hopping occurs. Since both the umklapp and 
the interchain hopping are relevant operators there is no 
weak coupling fixed point at the transition, and the RG 
(||) cannot be used to determine the critical properties 
of this transition. One can however get qualitatively the 
transition line as the position where the most relevant 
operator changes. Using (0) leads to 



(9) 



The resulting phase diagram is shown on Fig. |l| 

We now turn to the critical properties of the insulator- 
superfluid transition itself. Since it is difficult to extract 
from the RG flow, we can only analyse it in specific limits. 
The first one is the large transverse hopping t± ^S> t, U . 
On each rung there are two one-particle states, bond- 
ing or anti-bonding. The many bosons low lying states 
correspond to putting every boson in the bonding state. 
The chemical potential must be such as to ensure a 
particle- hole symmetry (around two bosons per "site") 
in order to have a constant density superfluid-insulator 
transition. This gives back the familiar problem of one 
species of bosons on a lattice with a commensurate fill- 
ing of two bosons per site, with an effective hopping t 
and an on site repulsive energy U/2. Since the effec- 
tive interaction is reduced compared to the single chain 
the large t± system can be superfluid even if the single 
chain is Mott insulating. The Mott-superfluid transition 
is a BKT one with an effective LU,narameter K e g = 2 
and a dynamical exponent z = lEloQ. To obtain the 
LL parameter K s of the original ladder one can com- 
pute the correlation functions. The one-particle green's 
function has a universal power law decay at the tran- 
sition (-^ ff (r)V>cff(0)) oc (l/r)^/ 4 ). Using @ leads to 
to the universal value K* — 1 at the transition. Since 
for fermionic ladders the large t± can lead to different 
phases than the small t± limit, it is important to check 
that an identical critical behavior is recovered in another 
strong coupling limit. We rewrite (^) as the bosonized 
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FIG. 1: Phase diagram as a function of i, t± and /x* 
at fixed U . /i* is defined as the distance in chemical po- 
tential from commensurability (i.e. one boson per site), 
that is /i minus a function of t and t j_ . This unimportant 
shift of /i is simply to keep straight axis in the figure. t\ 
is the critical value for a single chain. t2 is the critical 
value for an effective chain of 2 bosons/site and an in- 
teraction U/2 (see text). The wall of the critical surface 
correspond to a commensurate-incommensurate transi- 
tion with z = 2. The phase diagram as a function of K 
instead of t presents identical features since for U fixed, 
K and t are continuously related. 



form of a different microscopic lattice Hamiltonian than 
([j]), namely two coupled spin one-half chains: 

TT TjXXZ I TjXXZ JZ \ QZ QZ 



i 



H 



xxz 



Using the standard bozonized expressions for the spinst3 
one easily recovers (]?]) from ( |l0| ) with h x cx tj_, J\_ cx g u . 
J z , Jp are choosen such as to recover the correct bare 
Luttinger parameters K s and K a . The Hamiltonian fll(i| ) 
can be studied in the strong coupling limit h x ~ J*[ ^> 
>/a,/3 (i_L ~ U ^> t). Each rung corresponds to a four level 
system that may be easily diagonalized. Keeping only the 
two degenerate lowest levels we obtain a pseudo-spin one- 
half chain. If we denote its spin by S and its couplings 
by J the correspondance between the spin operators is: 



qz _ 

QZ _ 

Dp — 
Dp — 

cos(a) = 



cos(a)S' z 
S z 

— sin(a)/ 



(—1)" 1 * sm^a 

[(JiP + (2hfy 



(12) 
(13) 

(14) 

(15) 

(16) 



where / is the identity. Using these relations, (|10|) be- 
comes an XXZ Hamiltonian with the effective couplings: 



-J z 
J 



■h 



2hf 



J 2 



(17) 



A non polarized XXZ chain (corresponding to commensu- 
rate filling for the bosons) is in a LL phase for | J z \/J < 1 
and has a gap otherwise. For the pseudo-spin Hamilto- 
nian this means that as we increase the staggered field 
the model becomes gapless for a critical value of order 
h x oc JJ_. One recovers the BKT transition with an ef- 
fective universal critical LL exponent K c = 1/2 and a 
dynamical exponent z = 1. Using ([l2) and oc ipl^i 
for the bosons it is easy to check from (l|) that the origi- 
nal bosonic ladder has a universal LL parameter K* = 1 . 
The diagonalization on each rung shows that i has a 
non-zero mean value, which we interpret as an order in 
the field 6 a - These results for the transition are coher- 
ent with the previous large transverse hopping analysis, 
which gives some confidence that we indeed capture the 
correct critical behavior. The two above analysis extend 
easily to the situation where the chemical potential does 
not ensure particle hole symmetry In that case, this 
leads to either a magnetic field in the spin representa- 
tion or an extra chemical potential in the effective single 
bosonic chain limit. In each case the BKT iiiansition 
becomes a commensurate-incommensurate onellj, with a 
universal LL parameter K c s = 1/2 and a dynamical ex- 
ponent z = 2. The other properties at the transition 
(compressibility, etc.) both for the commensurate and 
incommensurate case can be obtained in a similar way 
than for the simple Mott transition for bosonscl. A sum- 
mary of the critical behavior is shown on Fig. 

Various physical observables can be computed. In the 
superfluid phase the symmetric modes are described by a 
LL with a parameter K s , whereas the 9 a mode is gapped. 
Thus most correlation functions involving the superfluid 
phase decay as power law. For example 



■;«-:c-k-„(i)')> x ( ~ 



(V4v!(r)^2(0)) cx ( ^ 



(18) 
(19) 



In the Mott phase <p s and cf> a are ordered so the corre- 
lations functions such as ([l8]) decay exponentially. One 
of the most important difference between the two phases 
are of course the transport properties. The Drude weight 
is zero in the Mott insulating phase, whereas it is given 
by T> = 2u s K s in the superfluid one, with a discontinuous 
jump at the transition. Far from the transition one can 
obtain K s from the RG. We use a two scale renormal- 
ization and cut the renormalization due to the umklapp 
when the the transverse hopping reaches a value of order 
one. This gives 



V = 2u a K s 1 



C_ 



-4K(2-K) 

L 4K -' -1 



(20) 
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ters could be extracted in a way similar to the one that 
was used for the single bosonic chainH. The gaps, super- 
fluid correlation functions and the Drude weight are prob- 
ably the most easily checkable quantities. For JJ arrays, 
as for the single chain the transition we predict should 
be visible as a metal-insulator transition in a transport 
experiment. One possible way to control t± could be to 
use a magnetic fieldEJ. 
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FIG. 2: The curve (a) is a schematic vue of the resis- 
tivity for the ladder in the mott-insulating state, with 
an activated behavior for temperatures lower than the 
charge gap A s . The curve (6) is for a superfluid ladder 
(shown for K < 1.5) (see text). The dashed part of the 
curves represents the cross-over region between low and 
high temperature. The insert shows the optical conduc- 
tivity (shown for K < 1.25). 



where C is a constant of order unity. Thus the Drude 
weight decreases as one approaches the transition by re- 
ducing t±. 

The temperature and frequency dependence of the con- 
ductivity can also be extracted from the RG. The high 
frequency behavior may be perturbativily computed from 
the umklapp scattering operator <7 u (cos2</>i + cos 202). 
This leads to cr(uO cx g 2 cu 2K ~ 5 , similar to the result 
for a single chainH. However for frequencies less than 
the gap in the antisymmetric mode, A a , the original 
umklapp operator has been renormalized nearly to zero. 
The dominant scattering operator is cos(v / 8</> s ) leading 
to (t(uj) oc gf l Lj 4:KB ~ 5 . In the Mott phase there is a gap 
A s in the symmetric mode below which these behaviors 
are cut, whereas it extends to u> = in the superfluid 
phase which has in addition the Drude weight at zero 
frequency. Similar behavior is obtained for the tempera- 
ture dependence of the conductivityQ: 
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p(T) cx T 2K ~ 3 
p{T) cx T 4K ~ 3 



T>A a 
T<A a 
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In the Mott phase one recovers the familiar activated ex- 
ponential behavior when T < A s . Both the dc and ac 
conductivity are shown on Fig. g. Note that for K s > 1.5 
the temperature dependence is non monotonous in the 
superfluid phase, whereas in the insulating one the resis- 
tivity would start to increase with decreasing tempera- 
ture even well above the Mott gap. 



The above predictions' could be checked either In nu- 
merical sirimkUiuris or hi (jxperimeriUl systems such as 

the J J arrays. In numerical simulations the LL parame- 
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